Directional (version 2.4)

MLE of the angular central Gaussian distribution: MLE of the angular central Gaussian distribution

Description

MLE of the angular central Gaussian distribution

Usage

acg(x)

Arguments

x

A matrix with directional data, i.e. unit vectors.

Value

A list including: A list including:

Details

There is a constraint on the estimated covariance matrix; its trace is equal to the number of variables. An iterative algorithm takes place and convergence is guaranteed.

References

Tyler D. E. (1987). Statistical analysis for the angular central Gaussian distribution on the sphere. Biometrika 74(3): 579-589.

See Also

rbingham, rfb, f.rbing, fb.saddle

Examples

m <- c(0, 0, 0, 0)
s <- cov(iris[, 1:4])
x <- MASS::mvrnorm(500, m, s)
x <- x / sqrt( Rfast::rowsums(x^2) )
mod <- acg(x)
mod
cov2cor(mod$cova)  ## estimated covariance matrix turned into a correlation matrix
cov2cor(s)  ## true covariance matrix turned into a correlation matrix